Search for: All records

Key Points A new semi‐analytical spin‐up (SASU) framework combines the default accelerated spin‐up method and matrix analytical algorithm SASU accelerates CLIM5 spin‐up by tens of times, becoming the fastest method to our knowledge SASU is applicable to most biogeochemical models and enables computationally costly study, for example, sensitivity analysis 

Abstract Soils store more carbon than other terrestrial ecosystems 1,2 . How soil organic carbon (SOC) forms and persists remains uncertain 1,3 , which makes it challenging to understand how it will respond to climatic change 3,4 . It has been suggested that soil microorganisms play an important role in SOC formation, preservation and loss 5–7 . Although microorganisms affect the accumulation and loss of soil organic matter through many pathways 4,6,8–11 , microbial carbon use efficiency (CUE) is an integrative metric that can capture the balance of these processes 12,13 . Although CUE has the potential to act as a predictor of variation in SOC storage, the role of CUE in SOC persistence remains unresolved 7,14,15 . Here we examine the relationship between CUE and the preservation of SOC, and interactions with climate, vegetation and edaphic properties, using a combination of global-scale datasets, a microbial-process explicit model, data assimilation, deep learning and meta-analysis. We find that CUE is at least four times as important as other evaluated factors, such as carbon input, decomposition or vertical transport, in determining SOC storage and its spatial variation across the globe. In addition, CUE shows a positive correlation with SOC content. Our findings point to microbial CUE as a major determinant of global SOC storage. Understanding the microbial processes underlying CUE and their environmental dependence may help the prediction of SOC feedback to a changing climate. 

Countries have long been making efforts by reducing greenhouse-gas emissions to mitigate climate change. In the agreements of the United Nations Framework Convention on Climate Change, involved countries have committed to reduction targets. However, carbon (C) sink and its involving processes by natural ecosystems remain difficult to quantify.

Using a transient traceability framework, we estimated country-level land C sink and its causing components by 2050 simulated by 12 Earth System Models involved in the Coupled Model Intercomparison Project Phase 5 (CMIP5) under RCP8.5.

The top 20 countries with highest C sink have the potential to sequester 62 Pg C in total, among which, Russia, Canada, USA, China, and Brazil sequester the most. This C sink consists of four components: production-driven change, turnover-driven change, change in instantaneous C storage potential, and interaction between production-driven change and turnover-driven change. The four components account for 49.5%, 28.1%, 14.5%, and 7.9% of the land C sink, respectively.

The model-based estimates highlight that land C sink potentially offsets a substantial proportion of greenhouse-gas emissions, especially for countries where net primary production (NPP) likely increases substantially and inherent residence time elongates.

 

The dynamics of soil phosphorus (P) control its bioavailability. Yet it remains a challenge to quantify soil P dynamics. Here we developed a soil P dynamics (SPD) model. We then assimilated eight data sets of 426‐day changes in Hedley P fractions into the SPD model, to quantify the dynamics of six major P pools in eight soil samples that are representative of a wide type of soils. The performance of our SPD model was better for labile P, secondary mineral P, and occluded P than for nonoccluded organic P (Po) and primary mineral P. All parameters describing soil P dynamics were approximately constrained by the data sets. The average turnover rates were labile P 0.040 g g⁻¹day⁻¹, nonoccluded Po 0.051 g g⁻¹day⁻¹, secondary mineral P 0.023 g g⁻¹day⁻¹, primary mineral P 0.00088 g g⁻¹day⁻¹, occluded Po 0.0066 g g⁻¹day⁻¹, and occluded inorganic P 0.0065 g g⁻¹day⁻¹, in the greenhouse environment studied. Labile P was transferred on average more to nonoccluded Po (transfer coefficient of 0.42) and secondary mineral P (0.38) than to plants (0.20). Soil pH and organic C concentration were the key soil properties regulating the competition for P between plants and soil secondary minerals. The turnover rate of labile P was positively correlated with that of nonoccluded Po and secondary mineral P. The pool size of labile P was most sensitive to its turnover rate. Overall, we suggest data assimilation can contribute significantly to an improved understanding of soil P dynamics.

 

Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

 

The terrestrial carbon (C) cycle has been commonly represented by a series of C balance equations to track C influxes into and effluxes out of individual pools in earth system models (ESMs). This representation matches our understanding of C cycle processes well but makes it difficult to track model behaviors. It is also computationally expensive, limiting the ability to conduct comprehensive parametric sensitivity analyses. To overcome these challenges, we have developed a matrix approach, which reorganizes the C balance equations in the originalESMinto one matrix equation without changing any modeled C cycle processes and mechanisms. We applied the matrix approach to the Community Land Model (CLM4.5) with vertically‐resolved biogeochemistry. The matrix equation exactly reproduces litter and soil organic carbon (SOC) dynamics of the standardCLM4.5 across different spatial‐temporal scales. The matrix approach enables effective diagnosis of system properties such as C residence time and attribution of global change impacts to relevant processes. We illustrated, for example, the impacts ofCO₂fertilization on litter andSOCdynamics can be easily decomposed into the relative contributions from C input, allocation of external C into different C pools, nitrogen regulation, altered soil environmental conditions, and vertical mixing along the soil profile. In addition, the matrix tool can accelerate model spin‐up, permit thorough parametric sensitivity tests, enable pool‐based data assimilation, and facilitate tracking and benchmarking of model behaviors. Overall, the matrix approach can make a broad range of future modeling activities more efficient and effective.